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In a paper published in this Journal recently, the analogy of an electric current 
in fluid motion was found to be what I there called the vortices of the vortices of 
the fluid. It soon became evident to me that the kind of motion there re- 
ferred to was a motion which was higher than vortex motion and yet had -some 
of its properties. 

The first portion of the present paper, treating of vector quantities in general, 
shows how to deduce the higher vectors from the lower, and also gives a method 
by which from three integrals any lower system can be deduced from any higher 
system. One of the most useful theorems, both for use in hydrodynamics and 
electro-magnetism, is that which shows how to replace any system of vectors by 
another system having the same external effect or the same internal effect, as 
the case may be. Another useful theorem shows how to deduce a system of 
vectors satisfying the equation of continuity from one which does not, though 
the principle of this theorem has long been recognized in the idea of cyclic con- 
stants. But I believe it will be found more useful as stated in this paper. 

Application is then made to the kinematics of fluid motion. 

Fluid motion must satisfy the equation of continuity, and so must be always 
cyclic. Otherwise there must be points somewhere in space where the fluid is 
created or destroyed. 

As Laplace's equation is a form of the equation of continuity, it is evident 
that this equation cannot, from a physical point of view, be satisfied at every point 
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of space, although there are some functions, like the spherical harmonics of posi- 
tive degree, that do so mathematically ; for these all give some fluid motion, 
often infinite in amount, at an infinite distance. So that such fluid motion can 
only be produced by the creation and destruction of fluid in regions infinitely 
distant. And so all fluid motion in an indefinite medium, without solid bodies, 
must be due, as it were, to the presence of vortex filaments or filaments of any 
of the higher motions* 

The n th order of motion always satisfies the equation of continuity, and hence 
the lines of motion are closed. All the infinite series of motions must exist some- 
where in space, but they may be confined to a point, a line, or a surface, and so, 
for all other points of space, only the n th and lower orders may exist. 

In this case, where the (w+ l) th order is zero throughout a certain region, 
the n th order will have a potential satisfying Laplace's equation throughout that 
region, though it may be cyclic. The motions of a fluid can be divided into two 
varieties, rotatory and translatory, the odd orders being rotatory and the even 
ones translatory. 

The motion of the zero order is ordinary irrotational fluid motion ; the first 
order is vortex motion ; the second order is what I have called the vortices of the 
vortices, or relative fluid motion, — in an element possessing this motion the centre 
moves faster than the sides ; the third order is one vortex filament inside another 
and rotating in the opposite direction : and so we can proceed. We thus get the 
idea that all the higher motions are, as it were, motions which take place inside 
the element as it drifts along. 

From this we see that when we deal with the higher motions, especially with any disconti- 
nuity in their distribution, we must use a smaller element than we otherwise should have to do. 

Thus, when we treat a single vortex filament, it is impossible to get the 
motion of the filament or the energy of the fluid without an integration across 
the section of the filament, however small this may be. This principle I con- 
sider of the utmost importance, and by recognition of it all that has hitherto 
puzzled us with respect to the motion of a vortex ring or the stability of fluid 
motion becomes a mere question of ordinary calculation. For, if we attempt to 
obtain the motion of translation of a circular vortex filament, we only obtain 
indeterminate results if we treat the filament as a whole. For I here show that 
the motion of the fluid exterior to the filament is the same for an infinite number 
of ways in which the vortex motion is distributed within the filament, and each 
of these ways has its own motion of translation and its own energy. 

A method is here given by which any distribution of motion of the n th order 
can be replaced by another distribution over a surface, so that the action on the 
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other side of the surface away from the original distribution shall be the same 
as before. This leads us, in the case of vortex motion, to the formation of hollow 
vortex rings, or solid ones of any finite cross section, and it is shown that one 
form of hollow vortex ring is permanent and probably stable. 

For I do not admit the validity of Thompson's proof that all discontinuous 
fluid motion is unstable. Illustrations are given of the distributions of motions 
of the different orders, and it is shown that any such series can represent a 
number of cases of fluid motion, seeing that we can take any one of the series 
of motions as the ordinary or irrotational motion. 

In treating of the action of forces on a fluid, we see that a single attracting 
particle can only produce fluid pressure, but no motion. Hence in an unlimited 
fluid forces having an acyclic potential can have no effect on the fluid motion, 
but can only produce a distribution of pressure. 

When a solid body is immersed in a fluid, the fluid is no longer continuous, 
but the surface acts like a diaphragm ; in this case we can have an acyclic 
velocity potential outside the body, but cyclic if the body was removed. But 
where the fluid is unlimited this cannot be so. 

The motion of a solid body in a fluid can always be represented by the proper 
case of motion in an unlimited medium, the condition being that the components 
of the fluid motion throughout a certain surface representing the surface of the 
body shall be the same as the components of the motion of an element of the 
surface of the body at that point, supposing the body to exist. 

As the classification of vectors is general, I here apply the method to forces, 
and so get the idea of forces of different orders. 

All forces which can produce fluid motion in an unlimited medium must be 
cyclic, and must contain all the infinite series of forces somewhere in space, 
though any given region may contain only forces of the » th and lower orders, the 
n th having a potential. 

Expressions are then given for the fluid pressure in terms of the motion of 
the fluid, one of which has already been given by Dr. Craig. 

I have not yet given any expressions for the energy of the fluid in terms of the 
higher orders of motion, though that in terms of vortex motion is well known. 

An immense number of applications of the ideas here set forth of course 
suggest themselves, but most of the results so far obtained are set forth in this 
paper, as I have not yet had time to descend to details. It is of course possible 
that some of the conclusions will not stand the test of time, but I hope that 
enough will remain to make the paper of value in all departments of physics 
where vector quantities are investigated. 
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The different portions of the paper have not been written in the order here 
given, and so the same idea is sometimes repeated in different parts. 

Neither is the notation very satisfactory, as different letters are sometimes 
used in the same sense, though in this case the connection is stated. I have 
changed the notation so often that I despair of perfection, and so publish the 
paper as it stands. 

General Theory of Vector Quantities. 

Let there be any distribution of a vector quantity, M , throughout space, 
and let its components be F ,G Q , and H . The dash over the letter signifies 
that these quantities do not satisfy the equation of continuity. Let us now 
subject these quantities to the following operations: — 

-n _ dHo dGo 

** ~ dy da' 

_ dF _ dff 
^ ~ da dx' 

_ dG __ dF 
M ^ dx dy ' 

dHy dGi 

Jfz ~'dy dz' 

^2 ~ dz dx' 

_ dGi _ dFj 
H *~ dx ~ dy' 

j-j dJT 2 dGi 

* s ~ dy da' 

etc. 

And let us continue substitutions of this nature indefinitely. 

We have thus derived from F , G , and lt an infinite number of new systems 
of vectors whose properties and relations to the original vectors are to be 
studied. In the first place, we notice that all these new quantities satisfy the 
equation of continuity except F , G , and If , and hence we have 

dF, dO, , dH, _ 

dx ~T~ dy "f" da U ' 

From these equations we may find, writing J for the quantity 

dF dG , dHo 
dx dy da 
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the following relations : — 

r * dx ^-^o? 

« = g-Alff p , 

and for the rest, 

Sn ~ ~~ A 2 ^T„_ 2 . 
Whence, writing xo for the quantity 

X ° = i fff'r (2 + if + S°) ^% & > 
we have, by Poisson's equation, 

^ = S JJJ T ^W^' - ^' 
^ = S /// T *%'&' - ^' 

^> = s J7X ~^ 2 ^w^ - w % 

and, in general, 

F n = -q J J J ~ dx'dy'dz', 

G " = kr J J J ~ J V ? dx'dy'dz', 

^ = £ fff ^T ****** 

From these, by differentiation, we have, after changing n + 1 to n, 

d- d- 

F n = br Jjj\ H " + l ij ~ G - + l ~ z \dx'dy'dz', 



when no Solid Bodies are Present. 231 

Let us now write 

F = F + #-"> 

m x n ' dx 

ri — n -i. ~x» . 

tr» — kr n T- dy > 



where 



v- 1 rrri\dF da dH n ) 



For all points where F n = G n = JBfj, = 0, — X B will be the potential of i^, 6r„, J2J,. 
From the form of these equations, F n , G n , and H n satisfy the equation of conti- 
nuity, while F n , G n , and H n do not. Hence every discontinuous system of vectors, 
F n , G n , H n , can be made continuous by the addition of the derivatives of a certain quantity 
having the value _ „ _ 

i rrri \.*f , do , dH) _ . _ 

4^ JJJ r \ dx + -dy + -di\dxdydz 

and the new system so formed satisfies the equation of continuity. 

To find a discontinuous system from a continuous one, we can draw surfaces 
which divide space into acyclic regions with respect to the kind of motion, and 
the discontinuous system will be distributed over these surfaces. 

If there is some necessary physical condition by which the vectors must satisfy 
the equation of continuity, then the discontinuous* system of vectors F n , G n , H n 
will evidently have the same effect in every way as the continuous system F n , 
G H , and H n , so that one system implies the existence of the other. 

This theorem is at the basis of all theorems of such a nature as those which 
express the energy of a magnetic system either by an integral throughout space, 
or by one throughout the magnets, or by one over the surface of the magnets. 

The value of ^„ can be expressed otherwise, as in Maxwell's Electricity, 
Art. 385, where it is shown that the integral can be divided into two integrals, 
one a surface integral and the other a volume integral. In comparing my 
formulae with Maxwell's, however, it must not be forgotten that my volume 
integrals are to be taken throughout the whole of space, whereas Maxwell's gen- 

* In this paper I nearly always use the term " continuous system of vectors " to indicate a system which 
satisfies the equation of continuity. 

But a distribution of vectors can evidently have a side discontinuity as well as an end discontinuity. The 
context will always indicate my meaning, especially as I use the dash over quantities which have end disconti- 
nuity and so do not satisfy the equation of continuity. 
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erally refer to the interior of some surface. Hence, in the case referred to, the , 
surface is at an infinite distance and the two volume integrals are equal. 
Consider the integral 

l rrri (dF dG dM H \ 
Xn = £ JJJ - r \~d^ + W + W \ dxd V dz 

taken without exception throughout space. 

This is equal to the sum of a surface and volume integral ; but if we take the 
surface where F n , G n , and H n are zero, the surface integral becomes zero, and leaves 

or we may write it 



i r r r m k cos <■ 7 . . 

Xn = ~ I? JJJ r 2 dxdydz, 



where M n is the resultant vector of the n th order, and e is the angle between it 
and r. 

If the point is within the region where F n , G n , II n exist, the integral will 
vanish at the lower limit, r = 0, and so the value can still be used within such 
space. For a region within which 

dF n . dG n . dH n 

dx dy dz 

the integral is much simplified. For it is only at the surface that this quantity 
has a value. Let the surface consist of an infinitely thin region -of thickness 
dv, within which M n decreases uniformly from its value inside to zero, then 













dK_ 


-_ dJ L? 


dv 














dx 


dv 


dx' 














dG«_ 


dG n 


dv 














dy 


dv 


dy 1 














dH n 


__dH» 


, dv 














dz 


dv 


dz' 




And 


since 


dxdydz = 


■ dvdS, 


we 


have 
















Xn 


i 

~ Ait 


■SS L 


cos 


dS, 



where <f> n is the potential of F n , Q n , H n , and is the angle between the radius 
vector, r, and the normal drawn outwards from the surface, and the inte- 
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gration is over the surface. Multiplying the equations by - dxdydz, and inte- 
grating, we see that we can write 

*. = h III %-' *w* ~ %" 2 ' 



*• = s //J ^ ^ w ~ 



dfn + i 

dz 



By differentiation and replacement of re + 1 by n, we have 

,i ,i 



F » = ^JjJ \ H ^+^- G *+^7u\ y ^\~^ — dry 

n 1 CCCSvi d \ f? d ^\i i i d(d£ n+ x dUi) 

G « = ^JJJ \ F n+i d j- Hn + 1 te\ dxdydz ~M~ dir ~~ d ^y 



dy 
where the value of £ m+a is 

&+« = -£ jjj ^ dxdydz 

or 4+2 = - ^?ffffffwr\^ + ^r + ^ 2 j**¥^ &%W- 

One of these integrations is to be taken throughout all space, supposing -F n +2> 
&n+2> f?n + 2 constant, and the other throughout space, supposing them variable. 
The first can be performed, for it is simply 

fffwr dxd y dz '> 

but we know that 

AV = -, 
whence we may write ** 

I I J jr dxdydz = — 2irr. 

These letters are connected by the relation that the R and r of the left-hand 
member and the r of the right-hand one form a triangle, the element dxdydz 
being at the intersection of the R and r of the left-hand member. The integral 
vanishes for the limit R = 0. 
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Hence £ l + 2 - §; J JJ rj-^j- + ~^r- + ~^r ^dxdydz. 

This can be put in the form 

& + a = — g; J J J -^» + 2 cos e dxdydz, 

where e is the angle between M n and the radius vector r. 

We can now develop the following general method of finding the lower 
vectors from the higher. Let us write 

O n = \ JjJF n rdx'dy'dz', 

P n = g J J J Gr n rdx'dy'dz', 

Q n = 2 J J J H n rdx'dy'dz'. 
Then we can write 

and the equations become 

•^»-2 4 T I ^ <A ^ «&. ^ <&, ^ dy T <fo J j ' 

r _IU 2 pj.iP.,^ , ^«\) 
tr»-2— 4^ j^^»^-<fy ^cte ^ ^ ^ <** / ) 

•On -2 ~ ±w I** v»^ dss \dx ^ dy ^ dz JY 

Let us now see what will happen if we multiply these by j^ dxdydz and 

integrate throughout space, R being a radius vector. 

The first terms become F n _ i} G n _±, and H n _ i , and we only have to consider 
the integrals of O n , P n , Q n . We have to obtain an integral of the form 

/ / / ■ j. dxdydz , dx'dy'dz'. 

Here we note that the radius vector r is the distance between the two points 
x ,y ,z and %', y', z', and so the integral reduces to 

ffff( x 'y' z '}fff \ *Z? dxdydz I dx'dy'dz'. 
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In this case, to obtain the integral with respect to r we make use of the fact 
that 

AV = n(ra + l)r"- 2 , 
and we thus find 

^ fffi dxd y dz = ~ (*+2)(»+3) rn+2 - 

The integral extends throughout space, and the r and R df the first member, and 
the r of the second member form a triangle, the point x ,y ,z being at their 
intersection. Hence, 

A 2 ffff(x',y',z')r n dx'dy'dz' = n(n + l) f fff(x',y',z')f l - 2 dx'dy f dz\ 

We have thus found F n _ i , G n _ i , H n _ i , and by repeating the operation s times, 
counting the original operation as one, we can find F n _ ie , (r„_ 2s , S n _ 2s , and 
F n _ 2s + i, Gr H -2s + i, H n _ 2s + i by differentiation of these. The result is 

°».» = lj^... 2 s JJJf^- 1 dxdydz, 
Pn, s = 1,2,3. ...2s J J J ®y s ~ 1 dxdydz, 

Qn,s = 17 2,3 2s J J J -S"^ 8-1 dxdydz. 

For the case s = 0, the coefficient of the integral becomes simply — 1 and £„ 
becomes — X«. We also observe that 

O nis _ t) = (AJ(O n>s ), 

Q n(s - t) = (AJ(Q niS ), 
and 

€n,s Air y dx ~*~ dy ' dz V 

In these the symbol (A 2 )' signifies that the operation A 2 is to be repeated t times. 
Hence we may write, in general, 
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and 



jp On,(s-l| 


di„,. 

dx ' 


-p 

ry ■*• n,(« — 1) 

tr»-2 S — 4t 


dL,. 
dy ' 


JJ ..._ Qn.it -1) 

■"» — 2s 4„. 


ds.' 


T7F _ dJl„—2s 

•^»-2s + l <fy 


dG n _ 2 . 

dz ' 


fcr»-2s + l— ^ 


dx ' 


-tin-2s + l— dx 


nF n -2. 
dy 



Thus we have expressed the whole series of vectors in terms of any one series, 
and have shown how to calculate the distribution of them all from the discontinu- 
ous distribution of any one order, and this by a single set of integrals, with the 
proper differentiation. 

The value of g n>s can be expressed simply by the integral 

&>• = S 1,2,1. ...2s jJJ M n^ s ~ 2 cos e dxdydz, 

where e is the angle between M n and r. For all space in which M„ is zero, 
F n ,G„, H„ have a potential found by making s = 0, whence 



L.o — — X» = ^ J J J n r T S ' <&tfy<fe. 



£„ satisfies Laplace's equation at all points where M n = 0, and this potential 
can be considered as due either to the distribution of M n within a surface, or of 
M n + 1 over the surface, the components of M n + 1 satisfying the equation of conti- 
nuity, but those of M n not. 

Indeed, we can consider any of the systems of vectors as the cause and all the 
others as the effect. 

When F n , G n , H n satisfy the equation of continuity within a surface, the 
integrals reduce to the surface integral of the next higher order of vectors. 

But in the above operation of finding F n _ 2s , G n _ 2s , S n -2s from F n , G n , H n , 
we have neglected the constants of integration, and 'we have now to consider 
their effect. 

We first observe that, as the higher orders are obtained from the lower by 
differentiating, they are perfectly determined by these equations. But as we 
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have to integrate in obtaining the lower orders from the higher, the result 
requires examination. Taking the equations for F n , G n , H n , and adding to 
them respectively v, ir ly ir n , we have 



F » = 4^ jjj ~T 1 dx'dy'M + 77, 
^» = 4^ J J J "i^ dx'dtfdz' + Wi , 
^ = K fff^T 3 - MUM + «*> 



Performing the operation A 2 on these, with respect to the point x', y', 2', and 
also noting that F„, G n , H n satisfy the equation of continuity, we shall have 

AV = 0, A 2 ir x = 0, AV U = 0, 

dir . (hri , d^u 

dx ~*~ dy * dz 

Now these quantities must not affect the values of F n + 1 , G n + 1 , H n+1 as obtained 
from F n , G n , H n , and so we must also have 

~dy dz ^' 

dir dir n 

~dz dx ^' 

dir 1 dir 

dx dy ' 

The solution of these equations is found in the values 

_dV _dV _dV 

n — dx ' "1 — dy' w n — dz ' 

where V must satisfy Laplace's equation throughout all space. And a term of 
this kind is the proper thing to add, no matter how many integrations are per- 
formed. But if this term satisfies Laplace's equation throughout all space, the 
cause of the potential V must be at an infinite distance. As this quantity is en- 
tirely arbitrary, and has no relation to the original system of vectors which we 
consider as the cause of the vectors of the lower order, we. may generally make 
it zero unless there is some other condition to be satisfied. But the integrals 
without any addition give us all the lower orders of vectors which are dependent 
on the higher. Hence no addition should be made to the integrals when we are 
searching for the effect of a cause. 
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Indeed, all such cases of a potential satisfying Laplace's equation throughout 
space can always be represented by a distribution over a sphere at an infinite 
distance, and in Physics it must always be looked upon in that way. In fluid 
motion the only quantities which we can add without destroying the original 
configuration is a general motion of translation and rotation of the fluid repre- 
sented by the equations 

7T = A + Bx + Cy + Dz, 
7ri = A x + B x x + C{y + DjZ, 
wii = Aa + B n x + C n y + D n z, 
where we have the relations 

B+C 1 + B u = 0; C U = A; B = B n ; C^B,. 

We can define V by saying that it is a solid harmonic of a positive degree, 
so that it becomes infinite at an infinite distance. 

We see, from these facts that if any quantities, F n , G n , H n , vanish, then all 
the others must have the values given by the above equations. 

Hence all the systems of vectors must exist somewhere in space, though it is 
evident that there may be regions where some of the orders vanish. When this 
is the case, and throughout the region the vectors of the n th order vanish, then 
all orders of vectors greater than the n th vanish also, throughout the same space, 
but not those less than this ; and the vectors of the (n — l) th order have a 
potential satisfying Laplace's equation. 

Furthermore, if within any surface we have a discontinuous distribution of 
the n th order of vectors, then we have seen that without the surface the vectors 
F n , G n , H n have a potential satisfying Laplace's equation, and consequently all 
vectors above the n th order will vanish outside that surface. 

If within or over a given surface the vectors F n , G n , H n form closed circuits, 
then the vectors F n , O n , H n , and of higher orders due to the above vectors will 
be zero without the surface, and the vectors of the (n — l) th order will have a 
potential. But in the surface, or without it, all the vectors exist. Conversely, 
if throughout any space the vectors of the n th order have a potential satisfying 
Laplace's equation, the vectors of the (n — l)' h order will be zero in that space. 

If within any surface we have a discontinuous distribution of the n tli order of 
vectors, it is evident that the external effect will be the same as the continuous 
distribution within the same surface of the (n — l) th or higher orders. Indeed, 
the effect will be the same throughout space. 
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We have seen that there are two methods of looking at a series of vectors, — 
one as if they were derived from those above, and the other from those below. 
Now these methods are similar to the direct and inverse methods of treating 
electrical problems. The inverse method invented by Green is by far the most 
useful in electricity, and it seems to me that a similar method will be useful in 
the treatment of vector quantities. Indeed, the method leads to most impor- 
tant results. Consider the equations 

_ dH n _ dG n 
-*» + ! ' dy dz ' 

_ dF n _ dH 

<r n + 1 ~ dz ~ dx> 

= dG^_dF 
■"» + 1 dx dy' 

The operations here indicate that we choose the vectors of the n th order, and 
find their cause in the vectors of the (n + l) th order ; whereas the direct opera- 
tion assumes the cause and finds the effect. This is similar to the case of elec- 
tricity where the cause of the potential V is discovered by the operation — j- A 2 

to be the density of electricity p, with its proper distribution. And just as 
we find many electrical distributions over surfaces which have the same ex- 
ternal effect, so we are able to find many distributions of the vectors F n + 1 , 
& n + i> H n + i, over surfaces which produce the same external values of F n , G n , 
and H n . 

Suppose that within a surface F n , G n , H n have a potential £ B)0 , and outside 
a potential ^ i0 , both of which satisfy Laplace's equation. To satisfy the equa- 
tion of continuity we must have at the surface, 

dv dv 

Our equations then show that F n + 1 , G n + 1 , H n + 1 are zero throughout all space 
except the surface. At the surface the equations become 

TP j ( d &.o _ dL,o \ dv_ _ f d£' nf0 __ d£ n _ \ dv 

-p« + i av — ^ dz dz J dy \ dy dy J dz' 

ri J — ( d &° _ dL,o \ dv_ __ ( d£' Hj) _ d£ n0 \ dv 

<jr n+1 av — ]^ dx dx J dz \dz dz J dx> 

TT rJ — ( d &° _ d Ahl\ * _ f d &.0 _ d Jn.o \ & 

Jln + iav — [^ dy dy J dx \dx dx J dy' 
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This distribution of F n + 1 , G n + 1 , H n + 1 will produce the potentials £ n>0 and 
^ )0 ,the one within and the other without the surface. The most interesting 
case of this discontinuity at a surface is found by making 

-*r = and *r = 

at the surface, in which case the surface is a stream surface for both f^ >0 and £ l0 , 
and the equation of continuity for the n th order is satisfied. 
Let us now make £ W|0 = constant, and we have 

7-, , <%»,o dv dj^p dv 

■Pn+iav— dz dy dy dz -> 
n a — ^"■° dv _ d &- dv 

tr n + 1 av— fa dz dz fa' 

77 A — ^ il _ d &-° d JL 

n n+1 av— dy dx dx dy -> 

which will apply to all surfaces as well as stream surfaces. 

If we had made £^ = constant instead of ^„ i0 , we should have obtained 
equations of the same form, but with the opposite sign. If, however, we say 
that the normal shall be drawn from the side of the surface containing the 
original source of potential to the other side, the sign will be the same. Thus 
our surface must be drawn so as to separate space into two parts, one of which 
contains the source, and our equations give us the values to distribute over the 
surface, so that, on the side opposite the source, the effect will be the same as 
from the original source. Hence, although we may express a proposition with 
respect to the outside and inside, yet we may always reverse the terms into inside 
and outside. This is exactly analogous to the proposition with regard to the 
distribution of electricity. If we distribute the vectors over the surface with a 
negative sign, they will just neutralize the effect of the original source on the 
side opposite to it, which is exactly analogous to the case of electrostatic in- 
duction. 

As to the total amount of the vector M n + 1 to be distributed over the surface, 
it is evident that the quantity in the two systems must be the same. - 

I here use the term quantity in the sense of surface integral across the section 
of the vector. Hence, if M' n + x is the original and M n + x the new value of the 
vector, we must have 

JfM' n+t dS' = JJM n+l dS, 
the integral to be taken over a surface cutting the vectors at right angles. 
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Having chosen our stream surface of the n th order, we now have to define the 
direction of M n + X on the surface. The equations express the fact that M n + l dv 
is in the direction of the intersection of the given surface with the equipotential 
surface, and has the value M n . The equipotential surfaces and two systems of 
stream surfaces can be found so as to. form an orthogonal system. But a stream 
surface in general is not necessarily one of them, but is only a surface containing 
the stream lines or lines of direction of M n . When we choose for our surface one 
of the stream surfaces of the orthogonal system, the intersection of the two 
stream surfaces is an equipotential line for M n + i, and the intersection of the 
stream surface and an equipotential surface is a stream line for M H + 1 . We have 
seen that the potential £' of the n th order can arise either from a distribution of 
vectors of the (n + l) th order in closed circuits or a discontinuous distribution of 
the n th order ; and that the (n + l) tb order was distributed in closed lines around 
the n th order. The relations are similar to those between magnetism and elec- 
tricity and between a magnetic shell and a current around its edge. In this case 
we can obtain a distribution of a shell of M n over the surface which is equiva- 
lent to the distribution of M n + 1 over the surface. The strength of the shell, 
M n dv, will evidently be found by adding all the shells together, and hence, if we 
draw any line I on the surface, 

M n dv = f M n + 1 dv cos 6 dl, 
but this is equivalent to 

M n dv — / M n dv cos 0dl — £ fl)0 + constant, 

because M n + 1 dv is equal to M n at the surface. 

Hence the strength of the shell is equal to the potential. Knowing M n + 1 dv 
originally, we can thus find M n dv, whence the potential is 



_ 1 f P M n d» cos* 



e being the angle between r and the normal to the surface. This integral must 
not only be taken over the surface which we have been considering, but also 
over diaphragms which we must draw so as to convert cyclic space into acyclic. 

Thus the theorem which I have just given merely shows us how to draw 
surfaces so that this integral taken over them shall be the same for all regions 
which can be reached without passing through a surface. 

It has already been known that the diaphragms could be drawn in any 
position, but this theorem applies also to the other positions of the surface 
boundary of the space under consideration. 
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We may evidently have two varieties of surfaces which divide the internal 
space into cyclic or acyclic regions. The first are obtained by forming, as it 
were, tubes around the closed circuits of the vector M n + 1 , which we can con- 
sider as the original source of ££ )0 . Distributing M n over this surface in the 
manner described, and removing the original source, the new system has the 
same external effect as the old, but makes f„ )0 constant inside. The second 
variety of surface is acyclic inside, and encloses the original distribution of M n + l 
as a whole. It is evident that this last variety of surface cannot be a stream 
surface for $' n>0 , unless there is some other source of M n outside the surface. Our 
process, then, merely replaces one of these sources for the points on the other 
side of the surface. 

Taking any cyclic value of the potential ^,o and drawing a tube of flow of 
any cross section, and applying this method to it, the distribution of M n + l over 
the surface will cause the potential inside it to become constant. Take away 
the external system now, and we shall have a potential equal to — ^ )0 inside 
the tube, and constant without. 

If we take any two of our surfaces and distribute + M n + l over the internal 
one, and — M n + 1 over the external one, according to the equations, then within 
the inside surface we shall have £ M)0 constant; between the two surfaces ^ has 
its former value ; and outside the outer one g n>0 is constant once more. In fluid 
motion this will give us the case of a mass of liquid revolving in a liquid at rest, 
with a core of liquid at rest. 

If we allow the original distribution to remain, and distribute — M n + 1 over 
the surface, then the potential inside remains the same and becomes constant 
outside. Thus the inside and outside of the surface are reciprocal, and we can 
always reverse the statement in this manner. 

Helmholtz has considered some cases of discontinuous fluid motion, and other 
writers have occupied themselves more or less with it. But most, if not all of 
them have considered the vortex surface between the two moving fluids as the 
effect^ whereas I here treat the vortices as the cause of the fluid motion, and 
have thus obtained a method of replacing the original vortices by one or more 
vortex sheets around them. This is the basis of the inverse method of treating 
hydrodynamics. 

To get a solid distribution of M n + 1 , we can place one surface within the other 
with its proper distribution so as to make what is required. If within a given 
surface any s} r stem of vectors are distributed so as to have a potential within that 
surface, then the . continuous system as obtained from this system will have a 
potential both within and without the surface, but the potential inside will not 
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be the same as that of the discontinuous system, but will be that minus ^ n0 . At 
the surface of discontinuity all the higher system of vectors will be collected. 

Let us suppose that we have the potential f B|0 within the surface, and ^ 
outside the surface, both of which satisfy Laplace's equation for the space in 
which they are used. We can get the conditions at any surface, but as they 
are complicated, let us take the surfaces, stream surfaces for both functions. 

Then at the surface 

dc Ki o ftgii.s p. 

dv ~ U > dv ~ U - 

According to our notation, £„ j0 is the potential of a higher order of motion than 
£' n<> , the first being the potential of F n , G n , H n and the last of F n _ 2s , G n _ 2s , 

■tin — is' 

If we should distribute — F n + 1 , — G n + 1 , — -S» + i, over the surface according 
to the previous investigation, we should have f B)0 = f n0 inside, and constant 
outside. If F n _ 2s + 1 , G n _ 2s + 1 , H n _ 2s + 1 , were distributed over the surface in 
the proper manner, Ave should have £' (>s = £ BS outside and constant inside. If 
the two existed together, there would be 

$n,o + constant inside, 
and 

£«,,s + constant outside. 

So that in this case we must distribute over the surface 



F n + l dv = - 


d£ n ,a dv c?f„, dv 
dz dy ' dy dz' 


&n + 1 dv = ~ 


d( n ,o dv «?4,o dv 
dx dz* dz dx' 




d£ n .o dv d£ nn dv 



and 



u n+1 av— dy dx -f- dx dy , 

w , <%„,, dv_ d£' n _ s dv 

J?n-2s + iav— dz dy >dy dz' 



d£', dv d£i„ dv 



*tW , 



G n -2s + idp— dx dz dz dx > 

it j — ^".» ^L d &-' — 

±i n _ 2s + 1 dv — dy dx dx dy -- 

Gauss's theorem applied to electricity gives us the amount of electricity in a 
surface, from the surface integral of the electric force. Let us now attempt to 
develop some similar method for this case. 
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Taking Poisson's equation, 

4 7rp = -A 2 F, 
Gauss's theorem gives 



i \ ~d~ d ® == — ^ w III pdxdydz, 



where the first integral extends over a surface, and the second over the space 
included within the surface. If we have a distribution of F n , G n , H n , within a 
given surface, then F n + 2 , G n + i , H n + i , are obtained from them by the equations 

F1+2 — — A 2 -F B , 

-"» + 2 = ' A. FL n . 
Hence we can write immediately 

J J -&? dS = - ffj F u + 2 dxdydz, 

Jj ~dv dS = ~~ JJJ G n + 2 dxdydz, 

Jf^dV dS = ~ fffffn + zdxdydz, 
As we have 

therefore we can also write 

fflHf dS = - fff'TJxdydz, 
where 

T _ dFn , d®n , dH* 
" n dx ' dy *~ dz 

To get a clear idea of the meaning of J n , I may remark that in the case of 
magnetism it gives us the distribution of the so called magnetic matter. In the 
present case the integral vanishes over any surface enclosing the original source, 
so that the algebraical sum of J n within the surface is zero, or the poles of the 
magnet are of equal strength. 

I have thus developed the general theory of vectors from a purely mathe- 
matical point of view, so that we may apply it to any subject involving these 
important quantities, and thus many excellent methods for use in electro- 
magnetism may be found. There are even some methods which give important 
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results in the theory of attraction and spherical harmonics, especially the theo- 
rem that 

A 2 ffff(x', y', O r n dx'dy'dtf = n (n + 1) ffff(%', y', *',) r n - 2 dx'dy'dz', 

A 2 being taken with respect to x ,y ,z, and not with respect to x', y', z'. This 
gives the potential with a force varying as r"~ 3 from that varying as r M_1 . 

But the most important and interesting application, is to hydrodynamics, 
and I shall then develop a few more points which might have been placed 
above. 

There is an important point connected with this theory which will be ex- 
tremely useful to us further on ; and that is that where M n exists and varies from 
one point to another, then we shall sometimes be obliged to use an element of 
the fluid in our calculations for the determination of M infinitely smaller than 
for M n . In other words, M will not be determined exactly within an element 
without going still lower and assuming something further within the element. 
Thus, in case we have an infinitely small tube having the vector M n in the direc- 
tion of its length, we find that there are an infinite number of ways in which the 
vector can be distributed within the tube and yet produce the same external 
effect. This may be unimportant in some cases, but a neglect of this principle 
will cattse many interesting problems in hydrodynamics to assume an indeter- 
minate form. Thus the problem of the velocity of a vortex ring needs to be 
calculated in this way. Assuming the distribution of finer vortex filaments 
within the substance of the vortex, and the problem becomes determinate. I 
shall thus show that the vortex ring moves according to that distribution. 

The Kinematics of Fluid Motion. 

The fluid is considered as unlimited and without solid bodies in it. The 
motion of a fluid is a vector quantity, and so the above theory must apply. 

Let jP , G , H , be the components of the fluid velocity,* and let the fluid be 
incompressible. Then we have 

dF dG dSp _ 
dx ' dy ~f" dz 

The vectors of the higher order, which I shall for simplicity call motions, must 
then be interpreted. The motion of the first order whose components are 

* We could evidently take any other order for the fluid velocity, and, in any problem, we can generally 
change the suffixes to the letters, and thus obtain new solutions. 



246 Rowland : On the Method of a Perfect Incompressible Fluid 

F x , G ly Hi, has been interpreted by Helmholtz to be vortex motion, and it 
remains to interpret the rest. 

The second order will be, as it were, the vortices of the vortices, but to in- 
terpret it better I shall proceed as follows : — 

Suppose that along the axis of x there is a line of motion of the n th order. 
Calculate from this the values of F , G , and H , and we shall have the true 
motion of the fluid. It will be found that in all motions of the even degrees the 
fluid moves parallel to the axis of x, but in motions of the odd degree the true 
motion of the fluid is in circles around the axis. Hence the motions of the even 
degrees are translatory motions and those of the odd degrees are vortical motions. 

Motion Of the zero order is ordinarily irrotational motion. 

Motion of the first order is the vortex motion of Helmholtz, and a line of 
motion of the first order is a vortex filament. 

Motion of the second order is what I have, in a previous paper, called the 
relative motion of the fluid, and a line of this motion has the fluid flowing forward 
along its axis with very great velocity and a less and less velocity as we proceed 
outward from the axis, until we reach a certain distance, where it is zero ; beyond 
this the velocity increases as we pass from the axis, and is infinite at an infinite 
distance in the impossible case of the line of motion being straight and infinitely 
long. In the possible cases where the line is closed, the motion of the fluid is 
finite everywhere, and is zero at infinity. 

Motion of the third order can be roughly conceived of as one vortex filament 
within another revolving in the opposite direction. 

Motion of the fourth order can be roughly conceived of as one line of motion 
of the second order within another larger and opposite one. 

And so, as we go upward, an element of the line of motion becomes more 
and more complicated. The motion of the fluid exterior to the line of motion of 
any order can be calculated by formulas as simple for high as for low orders of 
motion, though the integrations are generally somewhat more difficult as we go 
higher. 

Mathematically there are motions of the negative orders also, and they are 
useful in some calculations ; but the physical conception does not at first sight 
lead to such important results as the positive, though I hope to investigate this 
further. 

Let us now conceive of a closed line in space of any form whatever having a 
uniform distribution of one of the higher motions along it. 

As the curve is closed such a distribution will satisfy the equation of conti- 
nuity. 
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If the motion along the line is of the n th order, then there will be a distribu- 
tion throughout space of motion of the (n — l) th order and below, but no motion 
of the n th or higher orders. And the motion of the fluid in this case will be 
equivalent to that produced by a discontinuous but uniform distribution of motion 
of the (n — l) th order over any surface having its edge in the given line, the 
motion being perpendicular to the surface. Suppose now that the fluid in which 
this line is placed receives motion of the (n — l) th order, but no higher. In 
order that this may be possible, it must have a potential, and so will vanish for 
this region when we find the new values of the motion of the n th order for this 
region. But not so if we add motion of the n th order or higher. 

Hence we have the theorem that motion of the n th order is unaffected by 
motion of any lower order. 

We have seen that if within a surface there are closed circuits of the n th order 
of motion, then without the surface there is motion of the (n — l) tb order and 
below it, the (n — l) th order having a scalar potential. 

Hence, if two regions of. this nature be placed near each other, the motion of 
the n th and higher orders will be the same as before, but of the lower orders it 
will be changed- 

In considering the general relations of vectors, I have virtually shown that 
if the fluid has any motion, then the whole series of motions must exist sometvhere 
in space. But we have also seen that if one of the motions, say of the n th order, 
is confined to any given region, then all the higher motions will also be confined 
to that region, and the remainder of the fluid will have only the (n — l) tb order 
and below; and that the (n — l) th order outside a given region has a potential 
satisfying Laplace's equation, but the lower orders do not. 

Furthermore, I have also shown how to find the motions of any order outside 
the given region by integrating throughout the given region. 

Now suppose a given region has a distribution of motion of the n th order 
within it which is discontinuous at the surface. 

We have seen how to find the distribution of motion of the n th order without 
it, in order to satisfy the equation of continuity. We have seen that this distri- 
bution has a potential without the region, and that only motion of the n th order 
and below exists outside the region. 

Now let us suppose that the original distribution of the n th order within the 
region is also such that within that region a potential satisfying Laplace's equa- 
tion exists, but not the same as that without. Then it is evident that all the 
motion of the (n + l) th and higher orders exists only at the surface of separa- 
tion of the two regions, and that both within and without the region only the 
n lh and lower orders of motion exist. 
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Thus this surface constitutes a surface of discontinuity in the fluid motion of 
the n th order. 

Helmholtz has discussed the case of surfaces of discontinuity when the 
motion on the two sides was of the zero order, and thus the motion at the 
surfaces of the first and higher orders. 

We have treated this case very fully in the theory of vectors, and have there 
shown the exact distribution of the motions. 

The process deduced from this theory which is most useful in hydrodynamics 
is the following : — 

Let us have a distribution of the (n + l) th in closed curves so that it satisfies 
the equation of continuity. Then throughout the rest of space the motion of 
the (ft + l) th order will be zero, and the n th order will have a cyclic potential 
£„ , the source of the fluid motion being supposed to be the (n + l) th motion. 
Select any stream surface which will thus enclose some of the (n + l) th motion. 

Draw the equipotential surfaces for the n th order of motion so as to intersect 
the given stream surface in lines, and distribute over the surface, in the direction 
of these lines, motion of the (n + l) th order, so that 

M n + 1 dv = M a 

at every point of the surface. Then the new distribution of the (ft + l) th motion 
produces the same effect as the old. Therefore we can take away the old dis- 
tribution within the surface and replace it by fluid for which the potential of the 
n th order of motion is constant, and the external effect of the new system is the 
same as the old. 

By putting one stream surface within the other a solid distribution of the 
(ft + l) th order might be made. If we draw the equipotential surface in the 
manner that Maxwell has done, we can define the distribution of the (ft + l) th 
order as follows : — 




Fig. i. 



Let the figure represent a system of equipotential surfaces infinitely near 
together, intersected by two stream surfaces very near together. 
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Let d\ be the distance apart of two equipotential surfaces. Then we can 

write 

M n+1 dvd\ = M n dk. 
But we also have 

JV1 » dX ' 

therefore we have 

M n+1 dvdk = f^- dX. 

But from the method of drawing the figure the last member is constant, and 
consequently the first is. also. But the first is the total quantity of M n + l within 
the rectangle, therefore the proper distribution of M n + l is given by putting the 
same amount into each rectangle. Let us now attempt to get the total quantity 
of M n + 1 . By quantity of M n + 1 I mean the surface integral of its cross section. 
Thus, in this case, the quantity of M n + 1 is 



J M n + 1 dvd\ — — J -jY d\, 



the integrals to be taken around the section of the surface. The last integral is, 
in our notation, simply equal to the original quantity of (n + l) th motion. Hence, 
the quantity of (n + l) th motion distributed over the surface must be equal to the 
original quantity. The analogy of this to the case of electric distribution is 
apparent, the surface integral over the cross section of a vector taking the place 
of quantity of electricity. 

If we apply this to the case of vortex and ordinary motion, we shall have the 
case of liquid at rest enclosed within a vortex surface, with the motion of the 
exterior liquid irrotational and expressed by a cyclic potential. 

The vortices of the surface are constantly moving forward in the direction of 
the fluid motion immediately outside the surface.. 

If we wish a volume distribution of M n + 1 within the surface, we merely have, 
as I have shown above, to put surfaces within one another with the above rela- 
tive distribution, the strength of each surface being arbitrary. We have thus 
one other condition to fulfil arbitrarily. Let Q n be a function which is constant 
for the stream surfaces, or, in other words, the stream function. 

Then the solid distribution of M n +1 is represented by the formula 

Q 
M n + \ — dMll^(Qn\ 

where C is a very minute constant, and \Jj (Q n ) is a function of Q n . 

* The notation here used is that of the general theory given on p. 236, so that £ tt> o is the same as — X K used 
before. 
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The direction of M n + t is that of the intersection of the surfaces Q n = con- 
stant, and £, )0 = constant. 

Such a distribution of M n + 1 will produce a potential outside the surface 
which will be proportional to f re0 an d equal to it, if C has the proper value. 
But inside the surface the motion is dependent upon the form of »/» ( Q n ). 

The distribution of the vortices over all the surfaces need not always be in 
the same direction, but we may alternate as often as we please without altering 
the outer distribution of motion. 




Fig. 2. 

As an illustration let us take a circular vortex ring. The equipbtential 
surfaces and lines of flow for an infinitely thin vortex are represented by the 
preceding diagram, where the closed curves of somewhat circular form are the 
sections of the surfaces of flow, and the others are the sections of the equipoten- 
tial surfaces. 

The motion of zero order is distributed throughout space according to the 
cyclic potential $ n>0 , and the motion of the first order (vortex motion) is con- 
fined to the closed surface. 

If we should cause instantaneous vortical forces to act over a surface similar 
to the given one, the forces being proportional to the vortices in the first case, 
then the analysis of vector motions shows that these forces produce the same 
effect as the original distribution of vortical forces for all points outside, but no 
effect on the points inside the surface. Hence, the system formed by these forces 
will be the same as the system under discussion. The system thus formed will 
undergo changes which are investigated in the portion of this paper treating of 
the equilibrium of fluid motion. We have there seen that the whole system 
of changes can be treated from the drifting of the vortex filaments, the whole 
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dynamics of the problem being satisfied by making the product of the vortex- 
strength of each filament by its cross-section constant as the element moves 
along. In the present case, the elementary vortex filaments of the vortex sur- 
face are constantly drifting from one place to another, and they may thus finally 
obtain another distribution, in which case the form of the surface and its position 
may change. In the present case, if the fluid moves through the ring from left 
to right, and we refer to the upper section, the vortex strength tends to in- 
crease on the outer and right-hand side, and decrease on the inner and left-hand 
side. 

When there is any change of distribution, the ring as a whole will tend to 
revolve around the part which has increased the most, and so the whole ring 
tends, as a whole, to decrease in diameter and move forward, though the changes 
of form are not so obvious. 

As the ring moves, the weaker part constantly tends to go inside, and the 
stronger part outside, and so tends to a form of equilibrium, though this should 
be investigated more thoroughly. 

In treating a filament, we must descend to the elements of the filament itself 
and consider the motion of each secondary element as affected by the others. It 
is impossible to treat this filament or vortex surface as a whole. The fineness of 
filaments to which we must descend will depend upon the rate of variation of the 
vortex strength. Thus, when we are dealing with finite distributions of vortices 
without side discontinuity, the elementary filaments may be large, but when we 
treat filaments or vortex surfaces, the elementary filaments must be small in 
proportion. 

If we now consider the hollow vortex to be very small in cross section, com- 
pared with the radius, so that it is nearly circular, and so does not change its 
shape much, it is easy to see that the tendency of the elementary vortices to 
accumulate on the outside will tend to revolve the whole circle into a 'position 
farther on in the direction of the fluid motion inside the ring, and if we still 
further consider that there is a tendency to accumulate in the front part also, we 
see that the ring will tend to move forward and to contract in radius. If the 
cross section is very small, the form of the cross section becomes more and more 
stable. The radius of the ring finally becomes constant when the forward 
motion is such as to make the difference of velocity between the fluid in the 
inside of the hollow vortex and the outside, the same for every portion of the 
cross section of the surface, which is possible for a very small cross section, or 
for a large one if the shape is correct. 

In attempting to construct synthetically a hollow vortex which shall move 



252 Kowland : On the Motion of a Perfect Incompressible Fluid 

forward uniformly through space without change of form, we have simply to find 
a distribution of elementary vortex rings such that the fluid motion due to them 
combined with a uniform forward motion shall produce a stream surface along 
which the motion of the fluid is everywhere uniform as reckoned from a system 
of axes moving with the surface. 

Or, for the sake of ease in calculating, we can suppose the axes at rest and 
the fluid flowing past them. It would thus be perfectly possible to construct 
synthetically a hollow vortex ring of finite cross section and finite velocity of 
translation. 

In the case of a very thin, hollow vortex of circular cross section, the ve- 
locity can be obtained by finding the motion of the fluid in the interior due to 
all the vortex filaments of which the ring is formed. Or it is simply one half 
of the difference of velocity on the outside of the surface nearest and furthest 
from the axes of symmetry. 

The velocity of translation of very thin solid vortex rings can readily be 
found by calculating the velocity at the core due to the given distribution of 
vortex filaments constituting the ring. It will thus be found that for a given 
distribution of vortex filaments within the very thin vortex, the vortex must 
move perpendicular to the plane of its curvature and inversely as the radius of 
curvature, to insure stability. Thus, if we curve a circular vortex at any point 
outward, at first the curved part will move faster than the other part. Then 
each side of the bend will move outwards, and thus the ring will tend towards 
a circular form again, though it possibly oscillates around it. 

It is evident that whatever changes a hollow vortex may undergo the inner 
fluid can never mix with the outer, but the vortex surface forms a box through 
which not a drop of the fluid can escape however the surface may be twisted. 

Thomson and others have pronounced all such discontinuous fluid motion to 
be unstable, and to prove this he supposes a depression to be made in the surface, 
which he then supposes to increase indefinitely. That some surfaces are thus 
unstable no one can doubt. But I do not think that all are so. Besides, if we 
have a vortex surface drifting around in a fluid, the latter having only an irrota- 
tional motion, the depression which will form in the surface at any time will have 
to satisfy certain conditions which may neutralize the previous effect. I have 
given my reasons for regarding the hollow vortex of uniform strength of surface 
as the form toward which hollow vortices tend when distorted. But I am not 
yet prepared to give further results, but will wait to obtain an exact solution of 
the problem. It is very probable that, when disturbed, oscillations are set up. 
All the principles necessary for the solution are developed in this paper. 
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If it is found that the hollow vortex is unstable, it may still be possible to build 
up a solid stable vortex on the principles here set forth. The condition thus 
developed for a vortex surface which shall not change its form by the mutual 
action of its parts is that the strength of the vortex surface shall be constant at 
every part and the surface a stream surface. 

But in order that such a surface may exist, it must be remembered that the 
components of the vortex motion must satisfy the equation of continuity and the 
motion of the fluid inside the surface must be acyclic. This last condition is only 
necessary in order to avoid vortex motion in other parts than the sheet, which, 
however, may exist if we please. Having constructed such a sheet and calculated 
the fluid motion from it, it will have a potential everywhere without the sheet 
satisfying Laplace's equation, without vortex motion. If we construct a hollow 
vortex ring in this manner, with the condition that the surface be a stream sur- 
face, the fluid within the hollow part will have an acyclic motion, but such as not 
to affect the form of the surface. Thus the motion might be one of translation, 
and it would be thus possible to have a stable hollow vortex. The only form of 
surface which can exist without changing its form is thus symmetrical around the 
axis of motion, and has a cross-section which has not yet been investigated. 

In certain cases the changes through which a hollow vortex goes are periodi- 
cal, and it is a question whether, with a proper motion of translation, they are 
not always of this nature. 

Let there be a distribution F n , G n , H n inside a given surface. Then F n , G n , H n 
will have a potential satisfying Laplace's -equation at all points where F n = 0, 
G n = 0, and H n — 0, of the value £„ i0 = — ^»- 

or integrating by parts 

f-.o = Lffl W + ^ + %>»} ** - kfSS\\iz + f ' + w\*°*9**- 

If within the surface 

F * = lb' G n = Sjj> alld #» - ~te 

we have 

f».0 — krJJ^ C0S edS + k» 

where 6 is the angle between the normal to the surface and r, and the last term 
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disappears outside the surface where F„, G n , and H n do not exist. Whence out- 
side the surface 

£.0 = I //£ cos 0dS 
and inside 

? = hffti cos MS + <!>.. 

The integral is also evidently different outside and inside the surface. We 
evidently have also 

A 2 f«,o = outside the surface, 
and 

A 2 £«,o = inside the surface. 

Hence this is a case similar to the one before it, and the surface has a distribution 
over it of motion of the (n + l) th order given by the equations found there ; this 
motion is confined to the surface, and forms closed circuits on it in order to 
satisfy the equation of continuity. 

When motion of the (n + l) th order is so distributed over the surface that 
the component of the n* order of motion is zero in the direction of the normal 
to the surface, then, from what we have before proved, the motions of the 
(n + l) th order and higher orders will be confined to the surface, and will not 
appear in the remainder of the fluid. Thus in the motion resulting from a solid 
body moving in the fluid, the motion of the zero order is zero at the surface, 
and so there is no motion of the first or higher orders throughout the fluid, and 
the motion of the zero order has a potential. 

If the discontinuous system F n , G n , H n , whose component is M n , are dis- 
tributed over a surface with the resultant normal to the surface, dv being the 
thickness of the surface, and M n dv being constant over the surface, then we have 
seen that the shell so formed is equivalent to a line of (n + l) th motion around 
the edge. 

If dS is the cross section of the line, we must have * 

M n dv = M n + l dS. 

The linear distribution of M n + X is then equivalent to the surface distribution of 
M n . There will be a potential of the n th order throughout space, satisfying 
Laplace's equation : the motion of the (n + l) th and higher orders are confined 
to the boundary line of the surface, and are zero throughout the rest of the 
space. 



* This is proved in the portion of this paper relating to the energy of the fluid. 
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So that all space except the line surely contains no motion of the n th and 
lower orders, the n th only having a potential which satisfies Laplace's equation 
but is cyclic. 

The components of the (» + l) th order are readily obtained by differentiation 

of the equations 

M. +1 dS CI 

*—i k Jr as > 



Gn 



M n + l dS Cm _ 
_^ t±1 dS Cn 

-ttn-l— 4:TT J T aS > 

where I, m, and n are evidently the direction cosines of the element of the 
line ds. 

The scalar potential of the n th order is simply proportional to the solid angle 
subtended by the line at the point, and can be computed by known methods. 
The components of the n th order can also be simply obtained from the above 
components of the (n — l) th order by the proper differentiation. 

As the motions of the (n + l) th order and higher are confined to the line, the 
total motion of the fluid can be represented not only by a line integral of the 
(» + l) th order, but also by a line integral of any of the higher orders. The exact 
form I have not yet investigated. 

As an illustration of the methods here given, let us suppose that throughout 
a sphere of radius R there is a distribution of motion represented by F n , while 
G n and H n are zero. We shall then have 

O n , 2 — — 2i J J J r s dxdydz, 

Pn,i = 0, 

where the integral is to be taken throughout the interior of the sphere. We 
thus find 

Inside the Sphere. Outside the Sphere. 

^,2 = ^{35i2 6 + 105W + 21W-r 6 }; G£ >a = ^{3JR* + 422ZV+ 35r 4 }; 

•*»,2 ^Jm,2 -L »,2 %,2 "• 

Since we have A 2 r™ = n(n + l)r" -2 , we can write all the following quanti- 
ties from differentiating these : — 
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O n , 1 = ^ 1 1 5 1? + 1 tih* - r* J ; 



<l.i = 2j £r {*+**>]•> 



O B>0 = ^'{3i? 2 -r 2 }; 



P n,\ — Qn,\ — P »,1 — $n,l ; 



(K. 



A 2 0„, ( _ I) =4^„; 



3 r 

Pn,0 = H!»,0 == ■* n,0 ~ ^Jn,0 == U > 

a 2 o; (( _ 1) = o. 



From these we have 



-F> 



^-B 3 , 



^.^gfo^iP+UiZV-^cosfl; £.. = lio-{2iP + 5»*{coB0; 



£.,i = : §f{52P--»*}co80; 
£»,o = — "T rcosfl; 



__ ^E 3 cosfl 
£»,o — — 3 r 2 * 



Whence we can write the whole system of vectors as follows : — 



dL 



FJf 



*&., 



•^-4 = ^115^ + 10^^-^-^; ^-4 = #7{^ + 5 ^--ir; 



n d£n.2 

tm-4 — dy ' 



FT =- d ^- i . 
•"»-■* dz ' 



U„_4 — dy ' 

Jin-i— dz 



^„- 3 = 0; 



#„_3 — 



-* n 

30 



{5i? -»*}«; 



^-3 = 0; 

"■n-3 — 30r ) r 2 ~ °( Z ' 



^«-a = 3o"{5i? 2 -r 2 }y; 



*•..« — „. ««,,, 



'.-. = i , l^-»'l-Ti 



^»-2 — 



-H»-2 — 






-«»- 


-3 


— 30r \ i 


3 1" 


U 


K- 


-2 : 


F n E" 

~ Zr 


dx 


J 


G'n- 


-2 


dy 


'■> 




H' n . 


-2 




. 
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F 
tr„ - 1 — 3 2 ; 

TT F n 

-a»-i — — 3-y; 



""»— 1 — :- 



H'n-l ~ — 



3 r 8 ' 
F n R* y 



JP = W — ^-° = ? ^ • 

■fn J- n dx 3 •*» ' 

/-Y -- _ ^? — n . 

"» ~ dy ~ v > 

rr — — ^» — n . 



•*» «fe 3 r 5 ' 

n' — — ^k° = _ -^-g 3 3a;y _ 
"» rfy 3 r> ' 



7"T/ __^ "£tt,0 



F n R s 3xz 
3 r 5 ' 



The orders higher than these are only distributed over the surface, the next 
higher order being 

F n+1 dv=0, 

1-2 

G n + 1 dv = g F n £ > 



H n+1 dv=»F n 



B 



The next higher order is distributed, as it were, on each side of the surface whose 
thickness is dv. 

By giving to n different values up to n = 4, we can get various distributions 
of fluid motion, all of which satisfy the equation of continuity and are thus possi- 
ble forms of motion. When n = 4, the motion of the fluid can be considered as 
due to a uniform distribution of motion of the fourth order, F i} throughout the 
sphere, the motions of the lower- orders being distributed as given by the equa- 
tions. Or the motion can be considered as due to a distribution of F n + 1 , G n + l , 
H n + 1 over the surface of the sphere. When n = 0, the case is that of a sphere 
of liquid proceeding forward uniformly in the direction of the axis of x. 

If we replace the sphere by a solid sphere, the motion of the fluid outside will 
remain the same as before. 

The changes which a sphere moving in the manner described will undergo 
can be calculated from the drifting of the vortex sheet backwards, as shown 
further on. 

The sphere will thus tend to flatten in the direction of the axis of x and 
broaden out in the other direction so as to form a figure somewhat similar to a 
prolate ellipsoid. 
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It is to be noted that the whole system of motions which I have given above 
have only required one integration, and in the general case would only require 
three for the complete determination of all the vectors up to the 2 s order. 

As another illustration, take a distribution of F n along the axis of x, and make 
s = 2, as before. Let the area of the section of the small tube along the axis 
be a, and in the integration reject the part which becomes infinite. We 
thus find 



where 
whence 

and 



°n,2 = 64 9 l0 S 9- 



L,2 = 0, 

-T— = h A2 °-' = ¥ (9 s + 2 9 s ^g q), 
G n -i = 0> 

*'.-, = 0, 

£„_ 3 = -f- (1 + loggr)z, 

-^,-3 = — -f- (1 + logq)y; 



^»-2 = -^{3 + 21og 9 }, 

#n-2 = 0, 

#„-2 = 0; 
^.-.1 = 0, 
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The calculation of F^, Gf n , H n makes them zero for all points of space except 
the axis of x, just as it should. 

But we can also consider them as distributed over any surface enclosing the 
axis of x. 

Taking a stream surface for F n _ x , G n _ x , ff tt -i which will be any surface of 
revolution around the axis, the proper distribution of F n , G n , H n over the surface 
will cause the above values outside the surface, and a zero value of F n _ ly G n _ 1 , 
H n _ 1 inside. 

Thus, take a circular cylinder of radius 6. We must distribute F n a over the 
surface of thickness dv. Therefore we have 

2 irbF n dv = F n a, 

Beside this distribution of F n we can distribute it uniformly in the interior of the 
circular cylinder, or in any other way suggested by the equations. In the case 
of uniform distribution throughout the cylinder we can replace it by a surface 
distribution of F n+1 , G n+1 , H n + 1 and so on ad infinitum, the exterior distribution 
of velocities being the same, but the interior being different. 



The Action of Forces on Fluids. 

Let a system of forces whose components at the point x, y, z are X, Y, Z act 
on a fluid, and let p be the pressure of the fluid at the same point. The hydro- 
dynamical equations of Euler are then, 

1^- r_^_p^?_ n dF ° tt dF ° 
p dx~ -*° dt *o dx ^o dy M * dz ' 

1 ^E — v _ d Jh „ 2/t ^L° _ n d®o _ tt dGp 
p dy ~ *o dt *<> dx ^o dy ^o dz ' 

1^—7 _ ^9 w dH ° r< dH ° TT dH ° 
pdz~ Z <> dt -^o dx "^ dy M <> dz ' 

It is usual in treatises on hydrodynamics to consider cases where X , Y , Z 
have a potential, and it is there stated that vortex motion cannot be produced by 
such forces. But if we consider that forces which have a scalar potential are 
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such as are produced by direct attraction to or repulsion from points distributed 
throughout the fluid; we see that conservative forces acting in an unlimited me- 
dium can never produce any motion whatever, but only influence the pressure. 

Thus the equations as they stand contain much that is superfluous, and the 
motion will be the same in every respect if we differentiate in such a manner as 
to eliminate all portions of X, T, Z which depend on a scalar potential. 

Let us write 

then the equations can be put in the form* 

dW - dF,, _ 

~dx ~~ -^o "I" ~d7 — GqIIx— H Gi, 





dW 
dy 


— dt 1 

~ -*o "*" Ht ~ H* F x — F ff lf 




dW 
dz " 


— dS 

~ %o + ~dt = FqGx— GoFi. 


Let us now write 




_ dZ d¥ 
2L i~'dy dz' 

_ rfJSJ, dZ 
r i~ dz dx' 

dY 9 dX 
"l-~ dx dy* 



j , -rr dZ x dT x .„- d fdX,, ,dY ,dZ \ 

and also X a = -^ - -& = - A 2 J^+ Tx (^ + -^- +-% j , 

pr dXi dZi A2V" _1_ A. f^S _l_£Z? _l_^\ 

X *~ dz dx a Xo-T dy\dx ~T~ dy ~T~ dz J' 

_dj\ dX, - d /dZ, *Z d£\ 

** ~ dx dy ~ — ^ Z «> "*" dz\dx "•" rfy "•" dz J' 

etc. etc. 

Now I have shown in the theory of vector quantities that every system of 
discontinuous vectors can be replaced by another system which satisfies the 
equation of continuity, provided we can show a physical reason for the vectors 

• Given in a more restricted form by Lamb, in bis " Treatise on the Motion of Fluids," p. 241. 
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satisfying that equation. In the case under consideration, the pressure conducts 
the force applied at one point to another, and the whole system of forces so 
applied to the fluid must satisfy the equation of continuity. Hence for the 
system X n , Y n , Z n , can be substituted the system 

V XT' I Xn 

X n - Ji n + ^~, 

■p- __ fr _i_ *X» 



z„ - z. -t- ^ , 



where 



x " = -hJJJ\(^ + 'w+'^) dxd y dz 



When this substitution for the zero order is made in the original equations, the 
only portion of the pressure that will remain will be that which arises from the 
motion of the fluid and not from the applied forces. In this case we can simply 
write 

X 2 = -A 2 Zo; F 2 = -A 2 r ; Z 2 = -A%. 

Differentiating the third of our equations with respect to y and the second with 
respect to z,. and subtracting, we can write the first of the following series of 
equations ; the other two can be written from symmetry : — 

_ 8*1 dF dF dF ■ 

-* 1 ~ IF ~~ ±1 to ~ ^ ~dy ~ Hl ~di ' 

v S(ti _, dG ~ dG dG 

** ~ "&" — -* 1 &" ^ ~dj ~ Ml ~dJ' 

ry §H\ -j-, dJT dH dHg 

Zl "~ ~W ~~ •** ~to ~~ ***■ dy ~~ M * ~dz~> 

where the symbol 8 refers to the moving element and has the well-known value 

St — dt ^ ^° dx ^ ^o dp ^ n ° dz 
Performing the same operation on these, we have the equations 
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dF* 

'2 dx 



V °El IP "1 " _ /•' " - " XT 

-A-i — g£ -^2 J~ «"2 //« -"2 



r 2 % 



2 dz 



dFv dF\ 
dy dz 



dF\ dF ) 
dy dz 



- o \dG dG x dGjdG,, 
" r ^ + <fy rfz dy <fe 



+ 



rfiZs rfJTj 



r 2 = 



8G 2 
U 



dG 
'2 <& 



jFa ~1Z~ G-. 



dGo 
dy 



-tf 2 dz ~\- 4( 



dy dz 

dFt dF\ 
dz dx 



>■> 



dH x dH 
dy dz 

dl\ dFo] 
dz dx 



dG dG x dGi <W 
' dz dx dz dx 



SH 



dH 



Z * ~ U F* dx @* 



dy 



-^2-^ + 2 



dH a dH x 

' dz dx 

dF„ dF\ 
dx dy 



+ 

+ 



dG dGi 
dx dy 

dHa dHi 
dx dy 



dJIidHo 

dz dx 

dF\dFo) 
dx dy 

dG x dG,, 
dx dy 

dH x dffg 
dx dy 



When a fluid in motion is left to itself, the problem of the changes which it 
undergoes has never been satisfactorily solved, even in the case of a single 
vortex ring left in space by itself, and the principles to guide one in the solution 
of the problem have not been very satisfactorily given. 

When the forces acting on the fluid are zeroj we have 



St 
SG r 



dF e 



dF n 



-^1 dx + G l di, + -#1 



dy 



dFo 

dz '' 



St ~ *i dx "+" °"i dy + **i dz ' 
SH, „ dH„ . _ dH„ , „ dH ( , 



St — F i dx" + G i dy + ^i <& " 



dff. 



These equations contain the whole dynamics of the subject, and simply show that 
the product of the strength of the first order of motion by its cross section is con- 
stant. This principle then, applied in the proper way, contains the whole of the 
dynamics of a perfect fluid. 

The process indicated by the equations is as follows : Having given certain 
values of the components of the velocities, F , G , H , which satisfy the equation 
of continuity, we calculate from these the distribution of the first order of 
motion. 
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The variation of this first order of motion in each element as it drifts along is 
given by the above equations, or the variation of motion in any fixed element 
is given by these modified equations. 

The variation of the fluid motion can then be calculated from these. Again 
applying the method, we could get a still further change, and if we continued 
this step by step process indefinitely we might trace out the whole fluid motion. 
But we might originally obtain F , G , BT as functions of x, y, z, and t, so that 
for t = constant, they should satisfy the equation of continuity, and for t variable, 
the above equations. 

The result so found would give all the changes which a given system would 
undergo which was started at a given time in any one of the configurations. 
But it is to be particularly noted that in using these equations we must always descend to 
such small elements of the fluid that discontinuity is avoided. 

Thus we should never treat a single vortex filament by itself, but should 
descend to the still smaller filaments of which the vortex filament is composed, 
and apply our equations to these. 

In this manner the problem- of the motion of a single vortex filament becomes 
perfectly determinate instead of indeterminate as before. And I believe that a 
recognition of this subject will lead to extremely important results in this subject. 
It is, then, in the study of these differential equations that the final solution is to 
be obtained. 

We can always tell from the equations in which direction the system tends 
to change, and can thus base a theory of the stability of fluid motion on them. 
But I have not yet attempted to find any solutions. It is to be noted that the 
equations are satisfied at all points of space for which the first order of motion is 
zero, and so we can always confine our attention to points where it exists. 

But we are able to regard the matter from another point of view. When no 
external forces act on the fluid, the equations of Euler become 





1 dp _ 8F 
p dx ~ St 






1 dp _ SG 
p dy dt ' 






1 dp _ SH 
p dz dt 




m 


BG 


BR 


St' 


St' 


St 



where 



are the forces of acceleration of the fluid. Hence we have 



264 Rowland : On the Motion of a Perfect Incompressible Fluid 

dy St dz St ~ °> 

dz U dx U ' U ' 

d_ SG _ d_ BF __ 
<fcc &" tfy St 

These equations simply express the fact that in a fluid not acted upon by exter- 
nal forces the forces of acceleration are acyclic. One of the most interesting 
cases is that of a surface of discontinuity. In this case the equations assume 

the form 

/SW _ 8#o\ ± _ /8G£ _ 8GA dv_ __ 
\U U J dy \ U U J dz~ °' 

\U U ) dz \ St St J dx U ' 

(SGJ _ SGo\ dv _ /SK _ 8F \ dv _ 
\H~ St J dx \St St J dy 

To get the meaning of these equations let us transform them until the normal 
to the surface is in the direction of the axis of X. We then have 

S-ffi)' &H» » 

St St ~ U ' 

S£q 8£o p. 

St "St ~ U- 

And we also have for the continuity of the fluid 

K-F Q = 0. 

At all points of the fluid where there is no discontinuity, but the fluid velocities 
are obtained from a potential, these equations are satisfied. It is only where 
vortices exist that there is a possibility of the equations not holding. 

The equations show that there must be no discontinuity in the forces of 
acceleration at the surface and parallel to it, though there may be in the direc- 
tion normal to the surface. 

If we then have a surface of discontinuity of this nature, we must then add 
to the forces of acceleration in the interior of the surface others so as to make 
the system continuous at the surface in all directions except that perpendicular 
to the surface. This is the same thing as saying that forces must exist through- 
out the interior of the surface tending to change the configuration. 

These forces must evidently be acyclic within the surface, and have the 
proper value at the surface, and so are perfectly determined. 
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As an illustration of these methods of finding what dynamical changes will 
take place in a fluid, let us take the case of a single vortex filament along the 
axis of X. The fluid velocity will be inversely as the distance from the axis, or 

Go = F ia J> 
H = - F x a *, 

where q = Vy 2 + z 2 , and a is the sectional area of the filament. 

Let us take a stream surface bounded by x = b and x = and q = c and 
q = e, and distribute F x a over it according to our equations. Then, on taking 
away the original distribution of vortex motion, the fluid will move within the 
surface the same as before, but without will be at rest. We immediately see 
that the system is not in equilibrium, for the centrifugal force of the moving 
liquid remains unbalanced by the fluid pressure. We readily see that the whole 
ring will expand indefinitely. 

The second method expresses the fact as follows : For the surface q = c and 
q = e the vortex strength has of itself no tendency to vary. 

But on the plane surfaces we have 

^ = 

St U ' 

St 

SHi 

St 



°"1 q* "T" 


U > * 


^1 q* 


TT 2 _JH 

~ H l q* 



- F,a\a,^ - H,'-f\ 



Now the distribution of G x and 5i is in the radial direction, and so these equa- 
tions show that the direction of the vortex elements tends to change to one in 
a direction perpendicular to q. Thus we arrive at the same direction of change 
as before. 

By the other method of looking upon the problem, we see that forces of 
acceleration must exist in the stationary fluid, and so the stationary portion will 
tend to move. The direction is readily determined. 

We are now prepared to state what the conditions are that there shall be no 
tendency to change from its present configuration. 

For a surface of discontinuity the condition is simply that the surface be a 
stream surface and that the strength of the surface, Mdv, be constant, or that the 
fluid velocity on the two sides of the sheet be everywhere equal and opposite in 
direction. 
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This investigation has led us to look upon the subject of vortex motion 
from a broader point of view than before. For we have seen that if we reckon 
up the amount of this motion by the number of fluid elements multiplied by 
the strength of the motion in them, then the motion so reckoned is being 
constantly created and destroyed by the other fluid motions. But if we reckon 
the quantity of vortex motion by its surface integral taken across its cross 
section, or by its circulation, then the statement that it is indestructible by any 
motion lower than itself is perfectly gorrect. And this latter definition of the 
quantity of any vector is so important in the theory of the replacement of one 
system of vectors by another, that I propose that the term obtain general use. 

The equilibrium of fluid motion can evidently be either stable, unstable, or 
neutral. Thus an infinitely long cylinder of fluid revolving around the axis in a 
medium at rest would evidently be in unstable equilibrium. Thomson has 
given the criteria for such cases in terms of the energy of the system, and I am 
not yet prepared to discuss the subject further. 

But his conclusion as to the instability of all cases of discontinuous fluid 
motion I am not willing to admit. The case of the hollow vortex, it seems to 
me, shows that there can be stable forms. For the changes which the system 
undergoes are toward stability, seeing that the vortex distribution of the sheet 
tends to become uniform, which is a case of equilibrium. 

Again, the sheet can never be broken and the inside fluid mix with the out- 
side. Hence I am of the opinion that a hollow vortex ring is stable, and always 
tends to the form where the vortex strength of the sheet is uniform over the 
surface. But the subject should be carefully examined before a final decision 
can be reached. 

To determine the pressure in fluid motion we can use the equations 

1 dp_ __ SijJ 
p dx~ U ' 

1 dp _ _ W 
p dy ~ It ' 

1 dp _ _ 8J7 
P ds U ' 

whence we have 

A2 _ _ I <L 8 3. _i_ £ *&? i d ^ ) 

a P P\dx U ^ dy U W» St [' 

whence 

V ~^r JJJ r\dxW + dy--cT + d-z ~u\dxdydz. 
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This can be put in the form 

„ JL CfC\ dF <> d0 <> dF o dG * I dF a dH a dF a dff dGpdHo dG dff ) 
P~ 27rJJJ \dy dx dx dy ~^ de dx dx dz "+" dz dy dy dz ^Xdydz. 

We may also write a very interesting form from integration by parts 

p rrr\sF d} w d} w,*h 

P = ~I^JJJ \^dx- + -Udy-+^dz-\ dxd y dz - 

If the forces of acceleration have a potential throughout a given region, 
these reduce to surface integrals. If the forces of acceleration form closed cir- 
cuits, no fluid pressure exists ; that is, if the forces of acceleration satisfy the 
equation of continuity. 

In this case, since we have 



if we write 



d dF d^dGp d dff _ 
dx dt ' dy dt ' dz dt ' 

It o — *odx^ (x ody- i ~ M dz' 
„ _ j-, dG dG a dG 



we shall then have 
whence 



^~ J ■^'dx'^^'dJ- r ^'dz , 

rr> ijt dHp , n dlf , dH a 

1 ° — - t o~dx~ + (T »~dy- +I:l o~a^' 

v = £ Iff r | £ + £ + d i\ dxd y dz > 

P = ~rjff \tof« + M/^ + x/-£\dxdydz. 



dy 

If Mi is the resultant vortex motion and M the resultant ordinary motion, and 
we draw v in the direction of M , we can write the equation- in the form 

or, as we prefer to write it, 

A> = 2„j^-(f-) ! j. 

If points attracting as the square of the distance be distributed through space 
with a density 
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at every point, then the pressure would be the same as that due to the motion 
of the fluid. 

In this sense the above expression may be considered as the source of the 
fluid pressure. 

Let us now take the equations of page 260, and write 

1 = G Hi — H Gi, 
J= H F l — F H lf 
K= z F Gi — G Fi, 
whence we have, calling C a constant, 

p = G-% M ~ £ fffr J& + T y + & j dxd V dz > 
which can be put in the form 

p = C- \ Ml - £ fffl \M? -F F 2 -GG a - H Q H 2 ] dxdydz. 

We can also, from integration by parts, as the surface integral vanishes when 
the integral is taken throughout space, put 

p =C-%M$ + £fff\l^ + J^ + K d 4\dxdydz. 

This expression can be written in the form 

p=C-? 2 M2 + <l>, 

where <£ satisfies Laplace's equation at all points where no vortices exist In 
the expression as ordinarily given, (f> is the potential of the applied forces, but 
here it is given in terms of the vortex motion. 

Some of these expressions for the pressure are similar to those obtained by 
Mr. Craig, and published in the Journal of the Franklin Institute. 



